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We calculate for the first time the surface tension and 
curvature coefficient of a first order phase transition between 
two possible phases of cold nuclear matter, a normal nuclear 
matter phase in equilibrium with a kaon condensed phase, 
at densities a few times the saturation density. We find the 
surface tension is proportional to the difference in energy den- 
sity between the two phases squared. Furthermore, we show 
the consequences for the geometrical structures of the mixed 
phase region in a neutron star. 
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I. INTRODUCTION 

The possibility of different phase transitions taking 
place in the superdense interior of neutron stars has been 
the target of considerable interest during the last few 
decades, where pion and kaon condensation as well as 
quark deconfincment have been investigated. But only 
less than a decade ago was it realized that if the phase 
transition is of first order, then a geometrically structured 
extended region will form in the superdense interior of 
the neutron star, where the two phases are in equilib- 
rium 0j . The reason for this richness in structure is that 
a neutron star has two globally conserved charges, baryon 
number and electric charge, and two chemical potentials 
associated with these charges. Previous studies using the 
Maxwell construction could only ensure one chemical po- 
tential was common in the two phases, whereas the gen- 
eral phase equilibrium criteria by Gibbs 0] ensure ther- 
modynamical equilibrium for a system with any number 
of chemical potentials. The consequences are that the 
system is not locally charge neutral and a competition 
between Coulomb and surface energies are responsible 
for the geometrical structures. Moreover, the common 
pressure will vary with the proportion of the phases, and 
thus create an extended mixed phase region with struc- 
ture in the neutron star. 

For a first order deconfmement transition, studies of 
the detailed crystalline structure of the mixed phase re- 
gion have always been hindered by the lack of a single 
good model describing both phases. This is in contrast 



to the first order transition to a kaon condensed phase 
described in |||J , where both the normal nuclear matter 
phase and the kaon condensed phase are described by the 
same relativistic mean-field model, which allow us to cal- 
culate the profiles of all important quantities across the 
interface. From these profiles the Coulomb energy and 
the surface tension can be found, where, especially for the 
latter, only educated guesses were previously possible. A 
condensate consisting of negatively charged kaons is fa- 
vored in neutron stars because they, contrary to other 
kaon types, can replace electrons as neutralizing agents 

Detailed knowledge of the structure of a possible mixed 
phase region at densities above saturation is important, 
irrespective of which first order phase transition is re- 
sponsible for it, as it may have important consequences 
for transport and superffuid properties and rotation in 
the form of r-mode instabilities, non-canonical values 
of the braking index, and glitch phenomena in pulsars 

In the present paper we calculate the surface proper- 
ties, i.e., the surface tension and the curvature coefficient, 
for a semi-infinite slab of normal nuclear matter in phase 
equilibrium with a semi-infinite slab of kaon condensed 
matter and show the resulting crystalline structure in the 
central part of the neutron star. Since we assure com- 
pliance with Gibbs phase equilibrium criteria, the two 
phases cannot be separately charge neutral, though over- 
all the net charge vanishes. Thus for a infinite system we 
cannot explicitly take Coulomb interactions into account 
in the Lagrangian. However, it turns out that the typical 
radii of the geometrical structures are smaller than the 
Debye screening lengths of about 10 fm [fL5[ , and there- 
fore it is a reasonable first approximation to ignore this 
effect in the calculation of the surface tension. 

Section II contains a description of the non- uniform rel- 
ativistic mean- field model used to describe both phases. 
In Sec. Ill the surface properties are described, whereas 
the consequences for the crystalline structure in a neu- 
tron star are illustrated in Sec. IV. Finally our results are 
summarized in Sec. V. 
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II. NON-UNIFORM RELATIVISTIC 
MEAN-FIELD MODEL 

The relativistic mean- field model [16|,|l7j used here is 
described in detail in 0|. In this model the micleon and 
kaon interactions are treated on an equal footing which 
means both couple to the scalar meson a, the vector me- 
son and the isovector meson p, denoted <r, V^, and R p , 
respectively. The Lagrangian for the nucleons is given by 
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d^a - -m 2 a a 2 - U(a) - -V^V^ 



-mlRa ■ R», (1) 



where m* N = m N - g aN a, Vp V = d^Vv - d„Vp, Rp V = 
dpR u — d v Rfj,, and the scalar self- interactions are on the 
form U(a) — (l/3)bmN(g a jycr) 3 + (l/4)c(g a Arer) 4 , where 
b and c are constants, is the nucleon field operator 
and r)v is the isospin operator. 
The kaon Lagrangian is given by 



C K = V*K*V^K - rn§K*K, 



(2) 



where Vp = dp + ig uN Vp + ig pN f K ■ Rp and m* K = m K - 
g a K&- K denotes the K~ field and tk the kaon isospin 
operator. 

In the mean-field approximation only the O'th compo- 
nent of the vector fields and the isospin 3-component of 
the isovector field have finite mean values. The equation 
of motion for the kaon is 



{V^ + m£) K = 0. 



(3) 



In a semi-infinite system the kaon field is only transla- 
tional invariant in the transverse direction zj_ and can 
then be written as 



(4) 



where is the transverse momentum, lok is the in- 
medium kaon energy and 4>{z) is the kaon amplitude. 

For s-wave condensation, i.e. k = 0, the kaon field 
equation is given by 



V 



l 4> + (W + g U KVb + g pK Ro,3) 2 - m* K 2 ^j = 0. (5) 



In the semi- infinite case V 2 = d 2 /dz 2 , where z is the 
direction perpendicular to the surface. 
The kaon density is 

p K = 2 (uj k + gojK V + g pK Ro,3) K*K, (6) 

and only for an infinite system px = 2m* K K* K as seen 
from Eq. (||). 

The equations of motion for the position dependent 
meson fields become 
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m 2 <r 



gaN {-p s + brn N (g aN a) 2 + c(g aN a) 3 ) 
- 2g aK m* K K*K, (7) 



V 2 Vq - m 2 Vo = -g U N(Pn + Pp) + guKPK, (8) 



and 



V 2 i? ,3 - "7, 2 i? ,3 = -g P N(Pp - Pn) + g P KPK, (9) 

where the neutron and proton densities, p n and p p , and 
the scalar density p s are calculated in the Thomas-Fermi 
(local density) approximation. 
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Expressions for the local Fermi momenta ki are obtained 
from the Dirac equation for the nucleons 



Pi = gu>NV =f g P NRo,3 + ykf + 



*2 
l N > 



(12) 



where the upper sign is used for neutrons. The electron 
chemical potential p e = p n — p p . The kaon amplitude is 
zero unless the condition luk — Pk- — Pe is fulfilled. 

The model parameters gwlmi (i — a,u>,p), b, and c 
can be algebraically determined from the five bulk prop- 
erties of nuclear matter, which we take as: E/A = —16.3 



MeV, po = 0.153 far 
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= 32.5 MeV, K = 240 MeV, 



and m* /mN = 0.78. Because of the Laplacian terms in 
the meson field equations explicit values for the masses 
are required. We use for the u; and p their rest masses 
m u = 782 MeV and m p = 768 MeV, while the a mass 
is determined from the surface properties of symmetric 
nuclear matter as described in p7| , ^8[ , m a — 390 MeV. 
The kaon coupling constants for the vector mesons are 
determined from the quark and isospin counting rule, 



g^K = 5wAr/3 



and 



9pK — 9pN, 



(13) 



and the scalar coupling constant is fixed to the optical 
potential of the K~ at po from 



Uk(po) = -g<rK<j(po) - guKVo(po), 



(14) 
-120 



where we for the optical potential use Uk(po) z 
MeV. 

The Gibbs conditions for an infinite system of a kaon 
condensed phase and a normal nuclear matter phase to 
be in thermodynamical equilibrium at zero temperature 
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Pjv(Mn,Me) = -Pfcr(^n,Me), 
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(15) 



where Pjv and Pk are the pressures of the normal nuclear 
matter phase and the kaon condensed phase, respectively. 
These conditions combined with the global condition of 
electric charge conservation, 



qtotal = (1 ~ X)lN(Pn, Pe) + XQK^n, Pe), 



(16) 



where q denotes the charge of the corresponding phase, 
and the field equations (without the Laplacian terms) 
allow us to solve for the bulk values of the fields, densities, 
and chemical potentials in each phase for any volume 
fraction of kaon phase x between zero and one. 

In order to calculate the surface tension between the 
kaon condensed phase and the normal nuclear matter 
phase, the profiles of the fields and densities have to 
be determined across the interface between semi-infinite 
slabs of each phase. This is done by simultaneously solv- 
ing the four coupled differential equations for the K~ 
and the meson fields Eqs. (^,0]|) through a relaxation 
procedure, where initial guesses for the different profiles 
are relaxed to their equilibrium values. The boundary 
conditions at ±oo are provided by the bulk values of 
the kaon amplitude and meson fields for each phase at 
a fixed X- I n practice a 30 fm region with the interface 
placed approximately in the center is sufficient to fulfill 
the boundary conditions, see Fig. |l|. 



in the meson field equations, and only keeping the 
term in the kaon field equation. Differentiating twice 
brings the kaon equation to the form d 2 (j)/dz 2 = X/ 

The denominator 
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WJC + g^KVo + g p K Ro,3) — nb K 

vanishes in the bulk kaon condensed phase in contrast to 
the meson field equations, which take the form d 2 a/dz ~ 
Y/m 2 a (X and Y denote two complicated functions). 
Thus the Laplacian term is generally most important for 
the kaon field equation. In Ref. [|13|| the approximation 
was used to calculate the weak charge of droplets of kaon 
condensed phase to study effects on the neutrino opacity. 
We do both calculations and use the results from the ap- 
proximative procedure as initial guess for the solution to 
the full set of four differential equations. This procedure 
ensures rapid convergence. 



III. SURFACE PROPERTIES 

In the transition zone between two phases the pressure 
tensor is no longer isotropic as it is for the homogeneous 
bulk phases. For a plane interface with the normal in the 
z— direction, the pressure tensor can be split into a nor- 
mal and a tangential component. The normal component 
of the pressure P stays constant across the interface, cf. 
to Gibbs phase equilibrium conditions, whereas the tan- 
gential pressure component Pr(z) changes as function of 
z across the interface. The mechanical definition of the 
surface tension uq of a plane interface is |L9| 



to 




FIG. 1. Profiles of the meson fields <r, Vo, and i?o,3 and 
the neutron, proton, and kaon densities, p n , p p , and pk, re- 
spectively, across the interface. The boundary conditions at 
z = and 30 fm correspond to % = 0.107. 



It is a surprisingly good approximation and a sig- 
nificant simplification to neglect the Laplacian terms 
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(P - P T (z))dz, 



(17) 



where Pt(~°o) — Pt{°°) — P- 111 nuclear physics this 
definition of the surface tension is denoted Gibbs defini- 
tion 

Because we have assumed the leptonic species are ho- 
mogeneously distributed throughout the system, they do 
not contribute to the surface tension. The pressure of 
the hadronic species Pjj is given by 



Ph = p 



dp 



£H = PnPn + PpPp + PePk - £ff • 



(18) 



Therefore the surface tension of a plane interface between 
a kaon condensed phase and a normal nuclear matter 
phase can be written as 



dz(e H (z) - e ff:Ar - /i„ (p n (z) - p n>N ) 

- Pp (Pp( z ) ~ Pp,n) - PePk(z)) , (19) 



where the z dependent quantities at minus infinity take 
the bulk values of the kaon condensed phase and at plus 
infinity take the bulk values of the normal nuclear mat- 
ter phase. The energy density of the hadrons across the 
interface is given by 
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e H (z) = \ ((V,) 2 + mW) + \ ((W ) 2 + rnlV*) 
+ \ ((Vi? , 3 ) 2 + m 2 p R 2 Q 3 ) + 2 (uJK+g^KVo + g pK R^?4> 2 



7T 2 , 
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2 k 2 dk. 



(20) 



Notice, if the Laplacian terms are neglected except for 
the kaon field, the only consequence is that the squared 
gradient terms of Vq and i?o,3 change sign. 

The curvature coefficient is in Gibbs definition given 

by iB 



"(, = / dz(z - zo) (6h(z) - e H ,N - Mn {pn{z) - Pu,n) 

-00 

- fj, p (p p [z) - p p , N ) - ti e Pk{z)) , (21) 

where the surface location, i.e., the position of the equiv- 
alent sharp surface, 



zo 



I-00 Z P'( Z ) dz 



(22) 



LooP'( Z ) dz 1 

and the prime denotes the derivative with respect to z 



than the approximation. The curvature coefficient den- 
sity e c = (z — Zo)e s , and is therefore naturally more sen- 
sitive to the exact form, especially the width, of e s . This 
is also seen from the plot. The densities shown in Fig. || 
are for bulk parameters corresponding to x — 0.107. For 
X increasing from to 1 the width of e s increases from 
about 4 fm to 7 fm, while the peak value decreases. 
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FIG. 2. The surface tension and curvature coefficient den- 
sities, e s and e c , respectively, across the interface. The solid 
curves are for the approximative calculation and the dotted 
curves for the exact calculation. 

In Fig. @ the surface tension and curvature coefficient 
densities (i.e., the integrands of Eqs. ([n]) and (|2l|)), e s 
and e c , respectively, are plotted both for the full set 
of differential equations and for the approximation with 
only a differential equation for the kaon field. The differ- 
ences between the two surface tension density profiles are 
small. The exact calculation result in a slightly broader 
and more symmetric profile with a smaller peak value 



FIG. 3. The surface tension oq and curvature coefficient 
7g as a function of the volume fraction of kaon condensed 
phase X- The solid curves are for the approximative calcula- 
tion and the dotted curves for the exact calculation. 

The overall consequences for the surface tension and 
curvature coefficient are shown in Fig. ^, where <jg and 
— 7G have been plotted as a function of x- Notice, the 
curvature coefficient is negative. The surface tension 
decreases monotonically with increasing volume fraction 
of the kaon phase, whereas the curvature coefficient in- 
creases. The values of <jg from the exact and approxi- 
mative procedure deviate less than 0.4%. Thus the two 
curves lie on top of each other in this plot. This is in 
contrast to the curvature, where the difference is about 
33% between the exact and approximative calculation. 
Furthermore, the absolute value of the curvature coeffi- 
cient is comparable to the surface tension. This means 
the curvature energy is not negligible compared to the 
surface energy. But of course even higher order effects 
may mask the effect of the curvature energy as in the 
nuclear mass formula, where the apparent absence of a 
curvature term seem to be caused by higher order effects 
p3[ . For a spherical droplet of kaon phase with radius 
R = 5 fm in a system where x — > 0, the surface energy is 
reduced about 2^gI{&gR) = 18 %, whereas the surface 
energy increases, since the curvature is a signed quan- 
tity, about 40 % for a spherical bubble of normal nuclear 
matter phase in the x ~ > 1 limit. We have here assumed, 
that both ere and calculated from a plane interface 
remain unaltered. 
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In Fig. ^ the exact calculation of oq and —jg as a 
function of \ have been compared with the fits 

^K°- 00786 iiv7^) 2MeV/fm2 

/ _ \ 3/4 

*« - -°- 0976 (SvtS) MeV/fm ' (23) 

where the energy densities are in units of MeV/fm 3 . The 
fit deviate less than 5% for the curvature coefficient and 
less than 0.5% for the surface tension over the region 
where ex ~ £jv decreases from 682 MeV to 358 MeV, i.e., 
about 50%. It is not exactly trivial to understand why 
the fits, especially for the surface tension, are so good. 
For the surface tension, Eq. ( p3| ) is equivalent to saying 

(e(z) - e N )^ cx P N - P{z), (24) 

whenever e s = Pn — P(z) is not negligible. This is in- 
deed the case as long as cPe/dz 2 is negative but certainly 
not otherwise. We have not pursued the issue any fur- 
ther, but suspect this may be a generic property of phase 
equilibria described in mean-field theory in the Thomas- 
Fermi approximation. 
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FIG. 4. The surface tension ac and curvature coefficient 
7G, solid curves, compared with the fits given in Eq. (pj]), 
dashed curves, as a function of the volume fraction of kaon 
condensed phase X- 



Coulomb energies was always much smaller than the bulk 
energy. Furthermore, the Coulomb energy was, as usual, 
calculated under the assumption that the phase bound- 
ary is sharp. Figure 18 of that paper show the diam- 
eter D and the spacing S of the geometrical structures 
as a function of the radial coordinate for a star at the 
mass limit. As we have now shown, the surface tension is 
proportional to the difference in energy density squared, 
and it is between a factor of 2-5 smaller than what was 
guessed at in j3J. Knowledge of the density profiles of 
the charged particles across the plane interface allow us 
to calculate the Coulomb energy for a soft boundary. It 
turns out that the Coulomb energy is reduced less than 
15% for geometries with the typical small sizes encoun- 
tered here. This translate into about a 5% reduction of 
the sum of surface and Coulomb energies. Therefore it 
is a good approximation to treat the phase boundary as 
being sharp with respect to calculating the Coulomb en- 
ergy. Figure || is an updated version of Fig. 18 of Re f. [ff| . 
(Details of this type of calculation can be found in B4L) 
The differences between the two figures are minor. Over- 
all the sizes of the geometrical structures just decrease 
about 30%. This is due to the fact, that to first approxi- 
mation the location of the transition from one geometry 
to another is independent of <jg- Furthermore, the sum 
of surface and Coulomb energies scale only as crj/ 3 . 
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IV. CRYSTALLINE STRUCTURE IN NEUTRON 
STARS 

In Ref. Q the surface tension was assumed to be pro- 
portional to the difference in energy density between 
the two phases and the constant chosen rather arbitrar- 
ily, but in a way that ensured the sum of surface and 



FIG. 5. The diameter D (lower curve) of the geometries 
of the rarer phase immersed in the dominant phase and the 
spacing S between the geometries as a function of the radial 
coordinate r a . 

We have ignored the curvature term completely in 
these calculations, even though it is straight forward to 
include the curvature energy, see p5|| . However, if the 
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curvature energy is negative, as it is for kaon droplets, 
the sum of finite size and Coulomb energies will not have 
a local minimum but will decrease to minus infinity with 
decreasing size of the structure. This is of course unphys- 
ical and higher order effects will prevent it from happen- 
ing. These complicated higher order effects may, as previ- 
ously mentioned, even to some extent mask the presence 
of a curvature term. Moreover, the curvature coefficient 
is very sensitive to the exact profile of the phase bound- 
ary as we have seen, so its real value is uncertain for the 
rather small structures encountered here. Therefore we 
choose to consider the value of the curvature coefficient 
as some measure for the amount of uncertainty in the 
surface tension calculation, but not more than that. 



sition involving a kaon condensate. We do, however, re- 
alize there is still much room for improvement. 
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V. SUMMARY 



Little is known about the equation of state for nuclear 
matter above saturation density, but it is expected that 
at least a phase transition to deconfined quark matter 
happens at a few times the saturation density. However, 
this phase transition may not be the only one encoun- 
tered in neutron stars. A transition to a phase with a 
K~ condensate is also a possibility. If any of these phase 
transitions is of first order, a neutron star will have a 
mixed phase region in its dense interior, which is very 
likely to have some observable consequences. 

For the first time the surface properties in the interface 
between normal nuclear matter and kaon condensed mat- 
ter have been calculated, which makes it possible to study 
this mixed phase region in greater detail. Our calcula- 
tions are only a first approximation, there are a number 
of complicating aspects, which we have ignored or only 
treated approximately - e.g., explicit consideration of the 
Coulomb field which results in screening effects; the va- 
lidity of the assumption that the surface tension is the 
same for semi-infinite slabs as it is for small slabs, rods, 
and drops; and the importance of the curvature and even 
higher order terms compared to the surface term. Con- 
cerning the first two points, we can generally say that 
screening effects will reduce the Coulomb energy, likewise 
the surface energy will decrease due to the decrease in the 
surface tension when the system is squeezed. The reason 
for the latter is that the surface region will dominate a 
small system, so that only e s and not de s /dz is zero at the 
boundaries. For example a slab which is 10 fm thick and 
with a volume fraction of kaon condensed phase of about 
0.4, the surface tension is reduced about 12% compared 
to the surface tension of the infinite system (30 fm thick), 
while the absolute value of the curvature coefficient drops 
about a factor of three. These two effects pull in opposite 
directions with regard to the size of the geometry which 
minimize the sum of Coulomb and surface energies, thus 
the overall effect on the size is expected to be only minor. 

We have taken another step towards a better under- 
standing of the mixed phase region for a first order tran- 
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